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Abstract 

We show a new remarkable connection between the symmetric form of 
a quantum stochastic differential equation (QSDE) and the strong resolvent 
limit of Schrodinger equations in Fock space: the strong resolvent limit is 
unitary equivalent to QSDE in the adapted (or Ito) form, and the weak limit 
is unitary equivalent to the symmetric (or Stratonovich) form of QSDE. 

We prove that QSDE is unitary equivalent to a symmetric boundary value 

problem for the Schrodinger equation in Fock space. The boundary condition 

describes standard jumps of the phase and amplitude of components of Fock 

vectors belonging to the range of the resolvent. The corresponding Markov 

evolution equation (the Lindblad or Markov master equation) is derived from 

the boundary value problem for the Schrodinger equation. 
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1 Introduction 

The last two decades show a valuable progress in quantum probability theory and 
applications |ij-||. It was discovered that fundamental constructions of classical 
probability theory, such as central limit theorems, conditional expectations, mar- 
tingales, stopping times, the Markov property, and Markov evolution equations, 



* These paper is submitted to "Infinte Dimensional Analyses and Quantum Probability" and 
published partially in M, M, ||] 



have noncommutative generalizations U-|J. The quantum stochastic differential 
equation (QSDE) is a noncommutative generalization of the Ito stochastic equation 
suitable for describing irreversible Markov evolution in operator algebras 0— [11 



The QSDE and the Schrodinger equation describe a unitary evolution of a quan- 
tum system and its environment. The solution of the Schrodinger equation is a 
one-parameter unitary group U t , and the solution of the QSDE is a unitary co- 
cycle u(s,t), i.e. an interval-dependent family of unitary operators with one of 
two composition laws u(s,r)u(r,t) = u(s,t), s < r < t for the right cocycle, or 
u(t,r)u(r,s) = u(t,s) for the left cocycle. From a mathematical viewpoint, this 
structure difference between U t and u(s,t) is superficial and physically unobserval- 
ble, because u(s,t) is usually an interaction representation of some unitary group 
U t . s . 

A deep distinction between a QSDE and a typical Schrodinger equation is that 
the QSDE necessarily involves a singular component in its formal Hamiltonian, 
i.e. it can be regarded, up to unitary equivalence, as a Schrodinger equation with 
Hamiltonian operator perturbed by singular bilinear forms (see [0-[l |). 



Such Hamiltonians in Ti ® V (L 2 (1R)) appear as strong resolvent limits (r — lim 
or srs — lim) of self-adjoint operators H a = H <g> I + I cg> E + H^ depending on a 
scaling parameter a G (0, 1] such that 

H$ = iT ® A(f a ) +R®A + (f a )+K® A + (g a )A(g a ), 

where K = K* > 0, (R*)* = R, f,g e L 2 (R), f a (u) = /(aw), g a {u) = g{au), 
/(0) = g(0) = (27r) -1 ' 2 . The family of quadratic forms H a ^[h £g> ip(v)] — (h (8) 
i)(v), H a h®ip(v)) converges clearly to the quadratic form 

if*[/i <g> V(«)] = e M2 ((h,Hh) + \\h\\ 2 (^(v),E^(v)) 



+ (Rh, h)v(0) + (h, Rh)v(0) + (h, Hh) |v(0) |' 

with the singular component vanishing on the total subset consisting of vectors 
h®ip(v) eH(g)T s (L 2 (R)) such that v G C^(M), v(0) = 0, h G dom K n dom R n 
dom R* fl dom H. 

We prove, for commuting coefficients K, R, and H that, in the strong resol- 
vent sense, H a converges to a quadratic form with another regular and singular 
components: 



H4/Kg>^(v)] = e^ 11 [(h,H h) + \\h\\ 2 (il>(v),-Ety(v)) 



+ (h, Httyvip) + (h, H 2 h)v(0) + (h, H 3 h) \v(0)\' 



where 

H " = H ~ R 'lf^ R + R 'jf^ R ~ - G - Hl = R 'T^K - iL ' w ' 

Ill-defined operators A, A + , and A that correspond formally to singular quadratic 
forms 

A*[i>(v)} = e IH|2 £(0), At[i>(v)} = e IH|2 £(0), A^(v)} = e^ 2 \v{0)\ 2 

generate well-defined operator-valued measures 

M^T) = A{T) = f dtJ*AJ u M 2 {T) = A + {T) = f dtJ*A + J t , 
Jt jt 

M 3 (T)=A(T)= / dtJ?AJ t 

JT 

which, up to unitary equivalence, act as fundamental creation, annihilation and 
number processes in the Hudson-Parthasaraty framework: 
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A.(T)[rl>(v)] = e^l ^ dtv(t), At(T)[^(v)} = e" v " ^ dtv(t), 



A*(T)[i>(v)} = elHI 2 / dt\v(t)\ 2 , 



T 



where v(t) = .T-^^t^cj) is the Fourier transform of v(u). Set M (T) = mesT and 
denote by M(T) = J2o Lk <S> Mk(T) an operator- valued measure in the QSDE 

du(0,t) = u(0, t)M(dt+), M(dt+) = M(t,t + dt). 

One of our main results explains a connection between the unitary group Ut and the 
solution u(s,t) of QSDE: 

U t = e ltfi = s - lim e itSa = u(0, t)J t , (1.2) 

a— >0 

3 r — 

M(T) = Y^L k ® M k (T) = i / dr(J T Hj; - I <g> E) VT G #(iR), 

o •' T 

where L k = iH k (see equation (1.1)), and J t — I <g> e rfE , E = / duouoa + {uo) a{uo) 
is an environment energy operator generating an interaction representation of the 
unitary evolution 

u(s,t) = J s U t - s Jt- 



A remarkable observation is that the weak limit H = w — limH a ^ H = 
r — limi^Q, can be derived by a symmetrization of the QSDE: u(0, t)M(dt+) = 
u(0, t)N(dt), where M(dt + ) = M(t, t+dt) is the stochastic differential of the adapted 
equation, and N(dt) = N(t — dt/2,t + dt/2) is the stochastic differential of the corre- 
sponding symmetric equation. A relation between the symmetric differential N and 
the corresponding adapted differential M, proved in §2, has the form of an integral 
equation 

N(T) + I f M(dt+)N(dt+) = M(T) VT e B(M), (1.3) 

2 JT 

and an unexpected fact is that the weak limit H contributes to the measure N(T) 
of the symmetric differential: 

N(T) = i f dr(J T Hj; - I <g> E) or H = I®E+ . / J* T N(dr)J T . 

Jt i mes 1 jt 

Equation (1.3) implies the Hudson-Parthasaraty necessary condition for a solution 
of QSDE to be unitary §: 

L = iH s - \l*L, L ± = -L*W, L 2 = L, L 3 = W - I 

where L and W are related to operators K and R as in (1.1), and H s = H — 
i/2L*(I -W)(I+ W)~ l L = H - R*K(A + K 2 )~ l R. We would like to remark that 
the solution of (1.3) differs from solutions of the equation M(dt) = exp{N(dt)} — I 
derived in [|T^] for the symmetric operator- valued measure iN(T). 

Very technical assumptions sufficient for the unitary property of solutions of 



QSDE were obtained in the last decade basically by perturbation methods M, (15 



|1"8| , but symmetric operators responsible for the unitary property of solutions of 
QSDE were not discovered. The difficulties are connected with the violation of 
the group property by solutions of QSDE (see PH), and with the violation of the 
symmetric property of the form-generators considered on the domain of Fock vectors 
with smooth components. More precisely (see J2IJ), the formal generator of the 
Schrodinger equation, which is unitary equivalent to QSDE, reads as a dissipative 
operator, perturbed by nonsymmetric singular (in the sense of [I2|— 1|14|) bilinear 
form. 

In this paper we consider a class of QSDE which appears as an interaction 
representation of the strong resolvent limit of the Schrodinger equations d t i/j t — 
iH a %p t in the Fock space parameterized by the scaling variable a. The basic property 
of the limit resolvent is the existence of standard jumps of amplitude and phase of 
Fock vector components belonging to 1Z: 

(N + l)- 1 (I®A(5+)-W®A(5-))y = (L®I)y, V^eTl (1.4) 



where W = {2+iK)/(2-iK), L = 2i/(2-iK) R,V = {^ , • • • , *»(wi, ■ ■ ■ , w„), . . •}, 

n 

A(5±)*„(r) = lira V * n+ i(n, . . . , r fe _i, e, r fc , . . . , r„), N^f n (u) = n* n (w), 

\l/ n (ri, . . . , r n ) is the Fourier transform T W ^, T of n-th component of Fock vector \1/. 

The support of discontinuities of vectors from R coincides with the support of 
singularities of the quadratic form associated to the formal generator of the limit 
unitary group. Our main observation is that the generator of group derived on TZ 
(i.e. in the set of Fock vectors with the standard discontinuity (1.4)) 

H = H ®I+I®E+iL*W® A(<J_), E = T T ^ f dTa + (r)a(r) id T ^l T , (1-5) 

Jm\{o} 

is a symmetric operator, that is ($, H^) = (H$, \l/) V\l/, $ 6 U; its symmetric 
property can be verified independently under assumptions less restrictive than used 
for an explicit construction of the resolvent. 

The algebraic computations permit to derive the Markov master equation di- 
rectly from the boundary problem for the Schrodinger equation. 

The structure of the paper is the following. First, in §2 an algebraic analog 
of the Markov property is introduced for operator-valued processes in the Fock 
space. The processes with this property are called interval adapted. It is proved 
prove that the assumption (1.3) is sufficient to transform an adapted QSDE to 
the symmetric form. In §3 we construct an explicit solution of the QSDE with 
commuting coefficients and prove (1.1)— (1.2) for this particular case. A class of 
explicitly solvable Schrodinger equations in the Fock space is analysed in §4. A direct 
computation shows that QSDE in the weak Hudson-Part hasaraty form coincides 
with the interaction representation of the limit Schrodinger equation. Next, in §4 
we consider the properties of the resolvent of the limit Schrodinger equation, and 
prove that it ranges in a set of Fock vectors with components that satisfy a standard 
discontinuity conditions (1.4). In §5 the Markov master equation is derived from 
the boundary value problem for the limit Schrodinger equation. The main results 
of this paper were included in |2{J and |25| . 



This work, in its final stage, was inspired by intensive discussions with Prof. 
L. Accardi during a short stay of the author in V. Volterra Center in February 
1996; our interest in the problem of the parameterization of generators of QSDE by 
self-adjoint operators was initiated long ago by the papers |TJ], [JED], p|-|]r9| . 



2 Interval— adapted processes in Fock space 

Let H and He be Hilbert spaces. The Hilbert space of a quantum system and its 
environment has a form of the tensor product M = H®T S , where H describes states 
of the quantum system, and V s = V s (He) is the symmetric Fock space describing 
a state of the environment. We denote by (•,•), || ■ || a scalar product and a norm 
in the corresponding Hilbert space. The scalar product in V s is generated by the 
scalar product in He '■ 

(F,F') = F F' + J2-,(Fn,F^), F = {F ,F 1 ,...}eT s , F ed, F n e® n 1 H E . 

Let £ C V s be a total subset of coherent (exponential) vectors ip(f) = {1, /, / <S> 
/,...,/ <g> ... <8> /,...},/ E He (see |J, 0). A symmetric tensor product is defined 
on the total subset V of polinomial vectors Fj = f®i and F' n _^ = f'®' 1 ' 3 by the 
equation 

F®F'=\f -F* 0i F ®Fi+F^®Fi, ..., y Jts n (F j ®F' n _ j ), ...}, (2.1) 

L 3=0 ' 

where S n is a sum over all permutations of arguments in multipliers /, /'. This 
definition can be extended by continuity to the linear span of V and to SpanP = T s . 
Definition (2.1) implies the relation ip(f) ®ip(v) = ip(f + v) for coherent vectors. 
Hence, ip(v) = (SiiVK^T^) for any (finite) identity decomposition or complete system 
of projectors {tct-}t er C 13 (He), parameterized by subsets of a measurable space 
(T,T): 

n T = n T = n T , VT G T, n Tl n T2 = Tt T jt Tl = 0, VTi, T 2 : T\ fl T 2 = 0. 

The Fock space V s has a property of the tensor decomposition for any finite complete 
system of projectors. In particular, T S (C 2 (^R)) = ®jT T ., rf,. = T s (TC2(Tj)) for 
any finite disjoint partition {Tj}, and ttt = F* It J 7 for any unitary operator T, 
where T — M, He = >C 2 (-K); It stands for an multiplication operator by the 
indicator function of a measurable subset T. 

Later Mj- denotes a product H ® V s equipped by the system of orthogonal 
projectors {tit}t&t, and Mt is the subspace H ® T§,, where T§, = T s (titHe)- The 
operator-valued family u(-) : T — > B(IH) is called (tt,T)- adapted if for any T G T 
the following two assumptions hold: 

u(T)h <g> tjj(Tr T v) eH®T T , u(T)h ® ip(v) = tfj(Tr T av) ® {u(T)h ® ^(vf T w)}, 

where T a = M\T; that is the operators u(T) leave invariant vectors from Mt^ G if/^ 
(and the subset H? = H <S> Tf). In what follows the Fourier transform is used as 
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a unitary operator T , and the corresponding representation of multiplication by 
indicator functions play the role of projectors. In that case we call the family of 
operators {u(T)}, T e B(M) interval adapted (see [|16| ): 

{ip®i;{f), U {T)g®i;{v)) = eM(f,KT°v)}{<p®i>{K T f),u{T)g®iP{ir T v)), (2.2) 

with tp,g 6 H, 4(f), 4(v) & S, T a = M\T. Definition (2.2) means that the 
operator u(T) changes only the part of exponential vector from T^. Denote by 
B(IH T ) = B(7i) ® B{Yj) an algebra of all bounded operators acting as in (2.2). 
Equation (2.2) can be verified in the weak topology, and hence it holds for the 
adjoint operator u*(T) as a hereditary property of u(T). Relation (2.2) implies the 
commutativity of the components from rf a of a coherent vector and interval adapted 
operator families u(T), M(T) : 

{if j)(f), «(Ti) M{T 2 ) g 4(v)) = exp (/, % T »v) 

x {{u*^) cp ® 4(tc Ti /)} ® V(v?t 2 /), {M(T 2 )# ® ^t, v)} ® V>(% «)) , (2-3) 

where T a = M\ {UTj}, Ti fl T 2 = 0. An important example of interval adapted 
operators is creation, annihilation and conservation processes in the Fock space: 

(4(f),A(T)4(v)) = e^(I T ,v), (4(f),A+(T)4(v)) = e™(f,I T ). 

(4(f),A(T)4(v)) = (f,7r T v)exp{(f,v)}, (2.4) 

with 7T T = T$^I T (t)T^ t . 

Let a + (u>), a((j), u G M be creation and annihilation operator densities, E = 
/ uj a + {oj) a(u) duj e C(T ) (C(H) denotes a set of closed operators in a Hilbert space 
H) be a generator of the Journe unitary transform: J t = expjitE}. The group of 
unitary operators Jt4( v ) = ip(e tujt v) is uniquely defined, and from the canonical 
commutation relations we have: 

a+ (to) = J* a + (u) J t = e~ iu}t a + (u), J t a + {uj) J* = e lujt a + (uj). (2.5) 

Let W^JR) be the Sobolev space with the scalar product (f,v) = (f(u}),v(u))(l + 
co 2 ) p / 2 ), peR. Denote by A = A(l) = J M a(u) du, A + = A+(l) = J M a + (uj) du, and 

A = A + A, quadratic forms in W?(R), so that (4(f),A+Aip(v)) = 2vr f(Q)v(Q)e^' v \ 



(4(f), A^(v)) = V2^v(0)e (M , (4(f),A + 4(v)) = V^f(0)e^\ 

The quadratic forms A,A + ,A are not closable in r (£ 2 (iR)) and vanish on the 
dense set generated by coherent vectors {4(v) : v(0) = 0}. They are closable in 
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T S (W^{M)) (jg, 0) and define closed operators from T S (W}(R)) to T^W^iM)). 
Operators (2.4) can be expressed as operators, corresponding to densely defined 
closable quadratic forms: 

A(T) = -L f A t dt, A + {T) = -L f A+dt, A(T) = i- / ^<ft (2-6) 

V 27T ^T V27T JT Z7T JT 

with y4 4 = J t * A J t , A+ = J t * A + J t and A t = J* K J t . In the standard notation, the 
operator A(T) coincides with the annihilation operator A(It), where It(w) is the 
Fourier transform of the indicator function of a bounded set T G B(M). 
Let M(T) =J2 Lj <g> Mj(T) be an additive function so that 

M (T) = / ■ mesT, M^T) = A(T), M 2 (T) = A+{T), M 3 (T) = A(T), (2.7) 

Lj G C(TC), Mj(T) G C(r 5 ). Assume there exist a joint dense domain D G ?i for 
Lj and their products. Let g G -D , h &TC, \f\, \v\ < 1. Identities (2.4) generate the 
equation for correlators 

(/i ® V(/), M(r) g ® ^(v)) = e^ ■ £>, Ltf) • ((/)«*, J T u*) = O(mesT) (2.8) 

i=0 

with oti — 1 for z = 2, 3 and /% = 1 for i = 1, 3; «j = /3j = in other cases. It 
follows from here that the family M(T) is interval adapted in the sense of definition 
(2.2), and that the mean value is absolutely continuous with respect to the standard 
Lebesgue measure: 

%-lh ® Vlf), M(0, t) <? ® ^)) = e^ • £>, L ig )(f(t))<*v(t)* G L x (iR). (2.9) 
at i=0 

These observations are important for a rigorous definition of the QSDE in the weak 
form 0. 

Consider the operator-valued stochastic differential equation 

du(r,t)=u(T,t)M(dt+), r<t, s - lim u{T) = I (2.10) 

(see 0), where M(T) is an additive function on B(M) which ranges in the set 
of interval adapted closable operators wth a common dense domain T> C M. The 
strongly continuous interval adapted cocycle u(-) : B(M) — > B(H) is called a weak 
solution of equation (2.10) if for any ip G V there exist a limit 

/ u(s, r)M(dT + )(p = uj — lim J^ u(s,tj)M(tj,tj + i)(p, 

jeJiv(r) 



so that u(s,t)(p = f + Js u(s,r)M(dT + )ip, with tj = j2 N , Jn(T) = {j : tj G T}. 
From the bound (2.9) follows an existence of time derivatives of the cocycle u 
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u (h ® VCf ), u(r, t)g ® ^(«)) = £ (/» ® ^(/),^(r, t)L i9 ® ^(v)J /(£)^(£) ft (2.11) 

h E 7i, g G -Do, f, v € W^JR). Equation (2.11) is called a u>ea£; form of QSDE 
(2.10). In this section we aim to derive a symmetric form of the equation (2.10), 
that is to construct an operator- valued measure N(T) such that 



w 



lim e' 1 u{t) M(t, t + e)tp = w - lim e~ l u(t) N(t - e/2, t + e/2)ip V^ G V. 



We start with a suitable choice of operator- valued measures M(-) in equation (2.10). 
Let V be a total subset inU= Z/£, T = B(JR), n T = T^Jt^T^t- The main 
property of the projector ttt follows from the properties of the Fourier transform: 

e lwt Tt T e- iwt = n T+t , T w ^ t v= (2tt)- 1/2 / due^vtw), T e B(JR), t,uj G JR. 

Example 2.1. For equation (2.10) with bounded coefficients {Lj} we put T> = 
Ti ® Si, where Si C T s is a subset of coherent vectors if)(f) such that |/| < 1, and 
/ G C2(JR) is a finitely supported function. Si is total in V s (0) and X> is total in 
JH. It is convenient to put V — D <g) Si, for equations with unbounded coefficients, 
if there exist where is a dense joint domain D C Ti for all L^ and their products. 

Definition 2.1 27ie vector subspace of additive interval adapted functions M : 
B{JR) -► B(B^) «s /a&e/Zed 6y .A© ^/or any p e V 

(i) there exist an upper bound 

sup £|(^,M(7»|=/iT(^), (2-12) 

where T is a bounded Borelean subset of JR, {Tj} is a disjoint partition ofT, 
{ipj G JJJx a } is a uniformly bounded family of vectors from JH, sup- \\ipj\\ < 1, 
and ht is a finite measure such that \xt = fJ^rif) = if mesT = 0; 

(li) there exist a derivative 

= W,M(dt + )<p) = i f + 

at e^O 

which is a continuous function from JH — > C l ° c (JR) for any fixed <p G V. 
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Assumption (i) requires an existence of the weak absolute upper bound for integral 
sums for every natural N and partitions {Tj}; condition (ii) assumes an existence 
of the majorizing function (if), M(T)tp) which is absolutely continuous w.r.t. the 
standard Lebesgue measure. 

To fix a notation set ttt = J~t->wlTJ~u->t as a system of projectors, and B(M) as 
a Borelean a-algebra T. Denote 5n = 2 _JV = Ti, Tj = Tj ' — j 2~ N , Jn{T) = {j ■ 

Lemma 2.1 For every L G Av, V 9 ^ ^ an d the strongly continuous interval adapted 
bounded cocycle u(t) = u(0,t) there exist an interval adapted weak limits 



lim £ u(rf ) )L(rf\rfH) V = f u{r)L(dr + )^ (2.14) 

hm 2 £ «(f ) )i(^ ) ,r ) W + M^=/«W^ + M (2-15) 

jeJjv(T) 

called adapted stochastic integrals. 

Proof. We can apply the Lebesgue theorem on the dominated convergence to prove 
the equation (2.14) because there exist a uniform absolute bound (2.12) for sums: 

j n = E l(^X r i)£( r i> r j+i)v?)l <sup J2 l(^> L ( r i> r j+iMI <Mv). 

jGJiv(T) ^ j6Jjv(T) 

where tpj = u*(0, Tj) if), and the limit f T (ip, u*(r) <p) = lim e~ l (u*(r) ip, L(t, T + s)ip) 
belongs to L l ° c because of (2.13). Hence, lim J N = J T (ip,{p) = J T dr f T (u*(r) if),(p). 
Since f T (u*(r) •,</?) : M — ► L l ° c (M) is a continuous function for any (p E T>, there 
exist a unique element </?t G M so that (^, v?t) = ^t(V') V 9 )- Clear that ipx '■ B(M) — > 
iff is an additive function, and (■?/>, y?T n ) — >■ for any decreasing sequence {T n } such 
that fl„>iT n = because of (2.13). Therefore, ipx is a vector-valued measure, and 
we write ipx = J T u(t) L(dr + )ip. Thus, we prove (2.14). 
To prove (2.15) consider an identity 



' E 

iGJiv(T) j£J N (T) 



+ E «(,f>) (i(7f\ 7f> + W) " *0T + <W 7$)V 

iGJiv(T) V y 

As we noticed earlier, the first sum converges to integral (2.14). The second sum 
is absolutely bounded, and each term converges weakly to 0. Hence, the Lebesgue 
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theorem on dominated convergence implies a trivial limit for the second summand 

as N — > oo. n 

Let M^\ M( 2 \ ... G At>, for simplicity one can think that M^ k \T) = L k ®M k (T) 
(see (2.7)). 

Definition 2.2 Locally convex vector space Gv = Gv{M^\ M^ 2 \ . . .) C M denotes 
a closure o/SpanP C M, equipped by seminorms: 

Ptm = \m\+ sup (2.16) 

T^iik},^} lljm.es lj 

with 7\ n T 2 = 0, 4) G M T a, 1 1-0| | < 1, T a = M \ {T t U T 2 }, Tj e (0, £), te iR+. We 
assume that Gv is dense in M. 

The group J T leaves invariant the set Gv '■ Pt{g) — Pt(Jr9), Gv = J T Gv- 
Cocycle property of coefficients (in the sense of paper |F2TJ) is a simple conse- 
quence of the definition of J t '■ J*A(T)J t if}(v) = (lT,e luJt v)ip(v) = (lT+t,v)i(}(v) = 
A(T + t)ip{v), J*A + {T)J t ^{v) = A + (T + t)if?(v). This property of the coefficients 
M(-) of equation (2.10) generates a cocycle property of the solution u(T). Later we 
will assume that the measures M(-) G Av poses the same property: M(T + t) = 
Jl M(T) J t . It is clear that Gx> is dense in M under assumptions of Example 2.1. 

Theorem 2.1 Let the interval adapted cocycle u(t) be a bounded weak solution of 
(2.12), if G V and L, ML, M 2 G A v . If the semi-norm p T (L(0,r)(p) (2.16) is 
uniformly bounded for all r G M + , then the adapted and the symmetric stochastic 
integrals satisfy the equation 

J u{t) (L{dr+) + ^M(dT+)L(dT+)\<p = J u{r) L(dr)p (2.17) 

for every bounded set T G B(M). 

Proof. Let u(t) be a weak interval adapted solution of the equation (2.12). Consider 
a family of operators u s (t) = u(s)+u(s)M(s,t). The difference A u (s,t) = u(t)—u s (t) 
satisfies the equation A u (s,t) = f s f° u{t) M(dr + ) M(da + ) in the weak sense, and 
for every ip G M, ip G T> there exist a uniform bound 

\^,A u {s,t)L{r\ N) ,rfH) V )\ < ^-^sup ||tt*(r)^|| Pt-s(L(0, 2' N ) <p). (2.18) 

Since L(-) is an additive function, from the equation for u(t) we have 



u{tj) L(tj, r j+1 ) v = {u(t j+ i)L(tj, t j+1 ) - u(tj)M(tj,t. 



./• ' i+±) L ( T ii T j+± 
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-Ufa) M(t J} t j+ i) L(r j+ i,T j+1 ) - A u (tj, t j+ i)L(tj, T j+1 )jtp. (2.19) 

For every if) G M, tp G T> a prior bound follows from the definition (2.16) 

&jjf = K^^iTj) M(Tj,T j+ i) L(T j+ i,T j+1 ) (fi)\ <r\ SUp | \u*(r)lp\ | p n H, 

T6(0,t) 

and from (2.12), (2.16) we obtain 

&j% = \(ip, A tt (rj, r j+ i)M(rj, Tj + i)L(tj, r i+1 ) #)| 

<r x 2 sup Wu^rMp^mn))^, (2.20) 

T6(0,t) 

where the semi-norm p T1 (L(0, T\) tp) is uniformly bounded in N. Therefore, 

By £$? = <>■ 

Now it is sufficient to observe that for operators L, ML G Av the weak convergence 
in M take place by Lemma 2.1: 



lim ^ M(0,r i )L(r i ,r j+1 )v?= / u(0,r) L(dr+)</?, 
Bm X! «(0,r,)M(r„r J+ i)L(r„r ? . + i)v9 = -^ M (0,r)M(rfr + )L(rfr + ) V 9. 



X "Jjv(T) 



Hence, the weak limit is well-defined on T> by identity (2.19) and is called a sym- 
metric (in the Stratonovich sense) stochastic integral: 



u(T j+ i)L(Tj,r j+1 ) 



f u(r) L(rfr) = lim £ 

= j T u(t) (L(dr + ) + l -M{dr + )L{dr + )^ . (2.21) 

n 

Equation (2.21) is satisfied whenever the measure L(-) of the symmetric integral 
satisfies the equation 

L(T) + - f M(dt + )L(dt+) = M(T) VT G B(M) (2.22) 

2 JT 
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with the given operator-valued measure M(T) in the adapted integral. Assume that 

M(T) = ^£ i ®M i (T), L(T) = y ti l ®M l (T), 
o o 

M (T) = / • mesT, M X {T) = A(T), M 2 {T) = A + (T), M 3 {T) = A(T). 

In this case the equation (2.22) has the form J T {L(dt) + \M(dt)L(dt)} = M(T). 
Using the Ito multiplication table for stochastic differentials (see [[|) 





dA+ 


dA 


dA 


dt 


dA + 














dA 


dt 





dA 





dA 


dA + 





dA 





dt 















in the left part of the equation L(dt) + ^M(dt)L(dt) = ^ A 
system of linear equations with respect to £ k '■ 



Mi(T) we obtain a 



£ = £ + A4/2, & = *! + d£ 3 /2, C 2 = £ 2 + C 3 £ 2 /2, C 3 



+ C 3 £ 3 /2. 



Assume that the operator 2 + C 3 is invertible, and (2 + C 3 ) 1 : 7i — > dom/^. 



Then the solution of this system looks as follows: 
l 3 )-\ £ 2 = 2(2 + C 3 )- l C 2} £ 



2C 3 (2 + C; 



\-i 



2C 1 (2 + C 3 )-\ £ 2 = 2(2 + C 3 y 1 C 2 , £ = C - C x {2 + C 3 )- 1 C 2 - If we take 



C = iH + R* 



1 



2-iK 



R 



-G, £i = -L*W, C 2 = L, £ 3 = W-I, (2.23) 



as the coefficients {£k} satisfying the condition necessary for existence of a unitary 
solution of (2.1) (see @), then £ 3 = iK, £ x = iR* , £ 2 = «i2, 4 = %H, where 



K = 2i 



i+w> 



R 



2% 



i_,yyL are the coefficients of the symmetric operator H(T) : 

3 



iH{T) = Y,h®M k {T), M k (T) 

k=0 



T 



dt J* t M k Ji 



(2.24) 



with M = /, Mi = A, M 2 = A+, M 3 = A+A 

Since if is a self-adjoint operator, there exist a unitary operator W = (1 + 
i£f/2)(l - iiT/2)- 1 called a Cayley transform of if/2, and 2 + £ 3 = J + W = 2(1 - 
iK/2)~ 1 is an invertible operator. Furthermore, the operator C\ = iR*(I — iK/2)~ l 
is well-defined provided the operators K, R and R* have a joint dense domain D ; 
the operators C = iH + R*{2 - iK)~ l R) } C 2 = (I - iK/2)~ l iR = -WC\ are 
well-defined also. We call H, K, R, R* coefficients of the interval-adapted measure 
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M(T) = J T dt Yl £-k®JtMkJt) and the operators {—i£k} are called coefficients of the 
corresponding symmetric measure H(T) = —i j T dt J2 ik^JfMkJt- The assumptions 
on coefficients of the symmetric measure H(T) have a natural meaning: they are 
necessary for the measure H(T) to be a densely defined symmetric operator; and 
the equation (2.23) implies an implicit form of the same assumption. 

As we will see in the next Section, the relation between coefficients {£k} and {£&} 
can be also derived from the strong resolvent limit for a family of the Schrodinger 
Hamiltonians. Since the conclusions of §2 look unusual, we consider the weak and 
the resolvent limits for an explicitly solvable model in Fock space. By this reason 
we restrict ourself to the case of Hamiltonians with commuting coefficients. 

3 The strong resolvent limit of the Schrodinger 
Hamiltonians 

The weak and the resolvent limits of generators of strongly continuous unitary groups 
may be different in a quite elementary case. The next example gives a hint to our 
main result. Consider a family of unitary groups exp{itH a } = £7 t : 

Ul a) i)(x) = ip(x - t) vcpliX f drV a (x-t + r)\, x,XeR, ip e C 2 (M) 

where V a (x) = (27ro;)~ 1//2 exp{— x 2 /2a}, a G 1R+. Clearly V a (x) — > S(x) asa-> +0, 
and /q drV a (x — t + r) — > 7[ 0j t)(x), with an indicator function Ir(x) of the borelean 
set T. Therefore, the weak limit of the family of essentially self-adjoint operators 
H a = id x + XV a (x) generates a well-defined bilinear form on W^M) 

H, [<p, V>] = (if, H w ^) = i(<p, VO + A^(0)^(0) 

that corresponds to a formal operator H = id x + XS(x). 

On the other hand, the strong limit of the family of unitary groups U^ 

U t if>(x) = e lt ™ = lim Q Ut a) ^{x) = ^{x-t)e iXI ^ {x) = ^{x-t) Ue iX - l)I m {x) + 1 

implies the strong resolvent limit H which can be described by a formal generator 

H = r - lim# Q = id x + i(e iX - l)5(x) 
or by bilinear form 

H, [<p, ^\ = Jim j t (tp, U t il>) = i(<p, V') + i(e iX - 1)^(0)^(0). 
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Important observation is that H = r — limH a ^ H = w — \im.H a . The multiplier 
e tX — 1, an analog of the factor W — I = £3 in equation (2.23), appears in a manifest 
form in the equation for H. A similar phenomena happen in what follows. 

The range of the resolvent of the limit unitary group U t is a natural domain for 
the generator H. In this example it can be described explicitly by an equation: 

^00 /*oo 

Rrjjfx) = / dte-^Mx -t) + d(x)(e a - l) e -^ / dte _ "V(-t), 

Jo Jo 

where 8(x) is an indicator function of the half-line M + . The structure of the resol- 
vent shows that the functions from the domain of its generator H have phase jumps 
at the origin x = : lim x ^ + o R^ip(x) = e lX \im x ^^oR^ip(x). Hence the domain T> x 
of H consists of functions with a standard discontinuity at the origin: 

ip : i> G W%(m\{0}), lim tp(x) = e iX lim ip(x) 

and the operator H acts as id x if x 7^ 0. The left and the right limits exist at the origin 
for every function from V x because of the imbedding W^{IR \ {0}) C C(M \ {0}). 
Integration by parts and the identity (fi(x)ip(x)\+ = for functions from V x proves 

that the operator H is symmetric. The existence of a solution from T> x for the 
problem (H + ijj)ip{x) = f(x), x 7^ with the boundary condition as above for 
every / G £ 2 (1R), At = ±1 implies the self-adjoint property of H. 

Consider the Schrodinger equation dtipt = iH a %p t with the self-adjoint Hamilto- 
nian H a = H ® / + / CS> E + H^j depending on a scaling parameter a G (0, 1] : 

H$ = K® A+(g a )A(g a ) + R ® A + (f a ) + R* ® A(f a ). (3.1) 

Assume that coefficients H, K, R G C(7i) have a joint spectral family E x : 

H = J v x dE x , K = J A dE x , R = J p x e 1 ^ dE x , 

where u, p, $ are measurable real functions corresponding to operators H, K, and R; 
and functions f a {uj) = f(au), g a = g{a.uj) f, g are elements of the space L^~(M) = 

{/ : / G L 2 (M), f = Tzlt! e Lf(M)}, /(0) = g(0) = (2ir)-V 2 . We omit the 
indices a, when there is no conflicts. 

Denote by Pt = P t (A) a one-parameter unitary group in /^(-K) with the gener- 
ator N a (X) = uj + A|<? Q )(<7 Q | and prove that the solution of the Schrodinger equation 
U^ = exp{iH a t} acts as follows 

U{ a) h ® i){v) = J e lv ^ dE x h ®Jp t {X)v + ip x e^ f P s (X)f a ds) 
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xexp^e"*** j\f a: P s (\)v)ds-plJ* dr jT ds (f a , P r _ s (A)/ a )j. (3.2) 

For any bounded weakly measurable function P t : M — > B(Ti) the following equations 
are consequences of linear changes of variables: 



P r dr + 



P r4 _ f dr 



o 



r+t 



t+s 







P r dr, 



t rr 



JO 

t r r 



+ 



P r - P dr dp 

t+s rt+r ft ft+s 

i Jt 



t rs 



P r + P dr dp 



P r - P dr dp = / P r - P dr dp, 

JO Jt Jt JO Jt / ' Jo Jo 



t+s rr 



which are related to a group property of £/ t : the first equation generates the group 
property of the argument of the coherent vector, and the second one generates the 
same property of the normalizing factor in (3.2). 

The weak form of the Schrodinger equation can be verified by differentiation of 
the quadratic form <p t — (h® ip(v), U\ h ® i/j(v)) : 



<Pt 



J e l ^\dE x h,h)e^\^v,P t {\)v + ip x e^ J^ P s (X)f a ds ] j + 



* (/„, P s (X)v) ds - p\ f dr f (f a , P r ^(\)f a ) ds }, 
Jo Jo 



+ip\e 
We have clearly 

elHI 2 I (dE x h,h)(v x + \\(g a ,v)\ 2 



-W't 



t=Q 



+ 2Rep x e i ^(v,f a ) + j duj uj\v(uj)\ 2 \ = (h® ip(v),H a h® ip(v)). 

The one-parameter group C/ t is unitary because of basic properties of spectral 
operators E\ and the equation for the norm of an exponential vector: 



ip(P t v + ipe t9 J Psfds 



exp< \\v\\ 2 + p 2 



Psfds 



2 rt 

-2plme- 1 * / (f,P s v)ds 
Jo 



(3.3) 



5(a) 



Since P t is a unitary group, we have 



2Re / dr / ds P r 
Jo Jo 



drP s f dsP r 
o Jo 
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Hence, all terms in (3.3) besides |H| 2 are canceled out by the corresponding terms 

in the normalizing factor. Thus, we obtain 

Theorem 3.1 The unitary group of operators Ul = exp{itH a } acts as follows 

Ui a) h <g> V(u) = I e^dEx h ®il>(p} a \\)v + ip X f?* x f P<, (a) (A)/ a ds\ x 
x exp{ip A e-*** jf* (f a , P^(\)v) ds - p\ jf* dr jT (f a , P£i(A)/ a ) ds}. (3.4) 

Equation (3.4) contains four terms depending on a in the right hand side. The next 
Lemma "On Four Limits" (see \ffljj) justifies the limit of (3.4) as a — > 0. 



Lemma 3.1 Let P t (A) be one-parameter unitary group in L 2 (1R) generated by 
N a (X). Then for f,g£ L^~(M) such that /(0) = g(0) = -4= there exist limits: 

(1) ^ds(g a ,P^(X)f a )^(2-iX)- 1 ; 

(2) J 4 d8PM{\)f a (u) - e^l m (uj)(l - zA/2)- 1 ; 

(3) {g a ,Pt\X)v)^{l-i\/2)-^Z->tV\ 

(4) P t (a) (A) -► exp{iZ(A)7f (0)t )} = P t (A), exp{iZ(A)} = (2 + i\)/{2 - iX), 

with ttt a projector in L 2 (M) : 7?r = ^-k^tOO-P^-*-*- Limits (2), (3) are weak in 
L2(M), the limit (4) is strong in /^(-R). 

Proof. Limits (1-3) can be derived from the Taylor decomposition of the operator 
P t (A) in a power series J2 Cfc(«, £)A fc (as in |3|, [|) and sequential limiting procedure 
in each term as a — > 0. This operation maps the series into a geometrical progression 
that converges absolutely for |A| < 2. The Fatou-Privalov theorem enables to extend 
this result to the entire real line A G M, provided there exist a uniform absolute 
bound for the series in the upper half-plain ImA > 0. We can prove the required 
bound for the series corresponding to limit (4), but we do not have such bounds to 
prove limits (l)-(3). The method used below consists in application of properties of 
completely monotone functions (see [22|], Chap. XIII, §4) and the Duhamel equation 



for the group P t (A) : 

P t (a) (X) =e i " t + i\ I dse^-^g^ig^P^iX). (3.5) 
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(1) Set a t = a)- (X) = (g a , P} a (A)/ a ). As a consequence of (3.5), we have 

/ a s ds = ®t(ga,fa) + iX $ t -s(ga,ga)a s ds, (3.6) 

Jo Jo 

with $t(g, /) = J ds (g, e luJs f) = /i/, 9 (0, t), where /i/, 9 (-) is2 the probability measure 
such that 

^ g (T) = f ds f dtg(s-t)f(t), fx Lg (R) = 2irg(0)f(0) = 1. 
Jt Jm 

Denote by d p = d p a \\) the Laplace transform of the bounded continuous function 
a t . Since the equation (3.6) contains a convolution, we obtain an equivalent algebraic 
equation p~ l a p = ® p (g, f) + i\a p & p (g,g). The solution reads as 

a P = P$ P (g, /)(1 - i\p$ p (g,g))~\ 

where p§ p (g,f) = J °° e~ pt d(g,e luJt f) = J^f^ij)) is the Laplace transform of the 
probability measure fif >g . Since & t (g a , f a ) = jjf >g (0,t/a), then a p a \\) = d pa (\) and 

a Py A > ~ 1 „-\TT 73 - 1 , IM27T2 73 + ZA " 



1 - z A/i 9 , 9 (p) 1 + I A 1 2 fil g (p) 1 + I A 1 2 ii\ g (p) ' 

We recall that the class of the Laplace transforms of positive measures coincides with 
the class of completely positive functions. It contains the function (1 + |A| 2 p) _1 , and 
it is closed w.r.t. addition, multiplication and composition (see. ||22|| , Chap. XIII, 
§4). Hence, the real and imaginary parts of a pQ ,(A) are the Laplace transforms of 
bounded measures: 

gW(A) = f°° e- pt (A x (a- l dt)+i\B x (a- l dt) 

A (1R \ / i /,ff( iR +) 2 R (ID \ /*/,fl(^+K,ff(^ 



l + |A|X g (^+) 4 + A2' " +> 1 + |A|X 9 (^+) 4 + P 
It follows from here for any A € (F that 

/ a^\X) ds = A x (T/a) + i\B x (T/a) VT £ B(JR+) (3.7) 

Jt 

lim / a < f i) (\)ds = A x (lR + )+t\B x (lR + ) = (2-i\)- 1 . 

a^+0 Jo 

(2) Denote & a) (\,uj) = f* ds P s ^(X)f a (uj), and for v G L 2 (iR) set 

b[ a) (X,v) = f duv(u)ti a \\,u), <p t (v,f) = T ds(v,e^ s f). 



From equations (3.5) and (3.7) follows an identity 

b { t a) = ip t (v, f a ) +i\ ds p t - s (v, g a ) / dr (g a , P^ a) (\)f a ) 



= <Pt(v,f a )+i\J dstpt^iv^^^Axia-^t-s^ + iXBxia-^t-s))^. (3.8) 

The integral J T dte tuJt g a (uj) converges in L 2 (M) to e tuJt I( ^(uj) as a — ■> +0. Since 
L 2 (M) G L l ° c (M), the functions ipt(v, f a ) and (pt(v, g a ) converge to (v, e lu}t I( 0tt )). The 
function A\{a~H) + i\B\{a~H) is measurable in t, bounded uniformly in a G (0, 1] 
and converges at each point to (2 — iX)' 1 as a — > 0. Therefore, it is possible to pass 
to the limit in (3.8): ^ a) (A,w) -»• (1 + zA/(2 - zA))(u, e iw */( ,t)). Hence, for every 
v G L 2 (iR) 

lim f ds(v,P^(\)f a ) = (1 -iA/2)" 1 ^^^/^)), (3.9) 

a^O Jo 

that is w - km$ a) (\,tu) = (1 - zA^)- 1 ^/^*) (<*>)■ 
(3) Rewrite (3.9) in an equivalent form: 

lim f ds(v,P^(X)f a ) = -A- / ^^/ (0it) (T)dr. (3.10) 



■>(<*)( \\„.\ _ /■„. d(«) 



Since (g a , Pf (X)v) = (v , P_ t '(\)g a ), from (3.10) follows an equation 



lim jf' ds(g a ,P^(\)v) = -^-(F*^ v ,I m (T)) = J l-^(I m (r),^:^ T v). 

(3-11) 
The set of the indicator functions of bounded Borelean subsets in M is total in 
L 2 (M). Therefore, from (3.11) follows 

lim / dsr(t)(g a ,PW(\)v) = -*—(r(-),F^.v) 

ct^O J 2 — IA 

for every r,v G L 2 (1R), or w — lim (g a , P t (\)v ) = (1 — iX/2)~ 1 J 7 *_ > _ t v in L 2 (1R). 
From here we conclude that the norm of the weakly converging family 

g a ,t = (g a ,P t ia \\)v)eL 2 (lR) 

is uniformly bounded for every a G (0, 1]. 

(4) Note that the strong convergence of unitary operators follows from the weak 
convergence, and the weak convergence follows from the convergence ot correspond- 
ing quadratic forms. Hence, it suffices to prove that (v,P t (X)v) — > (v,P t (X)v) 
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Wv e L 2 (M). Denote c ( t a) = cf\\,v) = (v, P t (a \\)v). From (3.1) follows the iden- 
tity: 

C W = ( v , e ^) + *A f ds (v, e^-^g a ){g a , P^(\)v), (3.12) 



/o 

where the sequence g a>s = (g a ,P^(X)v) is uniformly bounded and weakly con- 
verging in L 2 (M). Set v(t) = T w ^>tV (w), an d v a (t) = (v,e luJt g a ). Since the function 
g{uj) = !F t _^Jg(t) can be decomposed into an integral, the square of the norm of 
$a — \\v — v a \\ can be represented as an integral of absolutely integrable functions 



$a = / . dtdrds (v(t — as) — v(t)) g(s) (v(t — ar) — v(t))g(r). 

Clearly, S a — > because it is possible to pass to the limit in this integral as a — ► for 
every v G /^(iR), g G L + ~(M). Therefore, we prove the norm convergence v a — » v. 

Now we can apply the equation (3) to the weakly converging sequence (g a , P^ a '(X)v) 
and pass to the limit in the integral (3.12): 

4 Q) - (y,e**{I+^£ dse—K-, s }v) = (v,e™{I + ^tfr m }v). 
In this way we obtain the weak and the strong limits 

lim P t {a \\) = 1 + j— -7f (0 ,t) = j— ^vf (0 ,t) + 7f(b,t) = exp{iZ(A)7T (0 , t )}, 



with the projector- valued family n? from B(L 2 (M)) , and I — n T = Ttx a - 
u 

Limits (l)-(4) describe what occur with the solution (3.3) as a — ► 0. Substituting 
(l)-(4) to (3.3) we obtain the unitary group U t = exp{iHt} = s — lim U[ . 



U t h®if>(v) = f e- Gxt dE x h^^U iZW ^°^e^ t v + ip^ x -——I {0 ,t)) 



x exp|ip A e-^^f^(7 ( o, t ),e^)|, (3.13) 

with Ga = — iv\ — p\/{^ — iX). The group property of Ut follows from (3.13) and 
from the properties of the Fourier transform: 

e 7r T e = 7r T+t , e J(o,t)(u>) = J( ,t)(u;J. 
Taking the time derivative 

1 ^7 

lim --{g <g> V(/), ^t ^ ® ^(u)) = H*[ 5 <g> V(/), /i ® ^(u)] 

t-*+o % at 
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where g, h & D C Ti, f,vE W^iJR) we obtain a bilinear form 

H*[9 ® W), h ® V(v)] = e (/ '^^ f (<?, F /i) w + (<?, if^)* u(0) + (g, H 2 h) H /(0) 



+(g, h) n J dujf(uj)g(uj) lu + (g, H 3 h) H /(0)v(0) ) , 
where 

■"0 = H — R - TT^R H" -R i TT^R = *"j •"! = -^ 



4 + if 2 4 + if 2 ' 2 -iff 

h * = iAk R > H *=o^r=Kl-W), W = l±^-. (3.14) 

2 — zii 2 — zii 2 — iK 

In what follows we assume that the operator G = — ii? + ^L*KL + |£*£ is a 
generator of one-paremeter contraction semigroup Wt = exp{— Gt} in TC such that 

D = dam H n dom L* L <Z dom G CdomL, G* <p + G<p = L* L<p \/<peD, 

H s = —H + ^L*KL is a operator symmetric on D, and D is core for G (see p3|). 

The formal operator expression describing the quadratic form H* reads as H = 
I® E + H (S)I + (2n)- 1 / 2 H 1 <g> A(l) + (2n)- 1 / 2 H 2 ® A+(l) + (2n)~ 1 H 3 <g> A+(1)A(1). 

Consider the function of a set u(s,t) = J s U t - s Jt- Since J^ip( v ) — ip(e~ 1 ' ujt v), we 
obtain from (3.13) 

u(T)h®i/j(v) = f e~ G * mesT dE x h ® ^(e iZ{X) * T v+ 

The m(T) family of operators is interval adapted, commutative for disjoint argu- 
ments, and satisfies the cocycle composition rule m(TiUT 2 ) = m(Ti)m(T 2 ), XinT 2 = 
0. The weak evolution equation for u(T) reads 

d(/i ® ij)(v), w(0, £)/i ® ^(f )) = z(/i ® ^(v), w(0, t)H(dt+)h ® ^(v)), 

where z#(r) = Af(T) = i/ r d£(J t HJ 4 * - E ® J) is well-defined operator-valued 
measure: 

H(T) = {H-R* I f r2 R+R* I f r2 R)^esT 

+ —^R®A + (T) + R*—?—®A{T) + i{I-W)®A{T). 
z — iK 2 — zA 

Hence, we obtain the following result. 
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Theorem 3.2 The family of solutions for the Schrodinger equation with Hamil- 
tonian (3.1) converges in L 2 (M) to the solution of stochastic differential equation 
(2.10) with coefficients (2.23): u(0,t) = s - lim a ^ UJ; a) J t * . Conditions (3.14) are 
necessary for the generator H of the limit unitary group U t to be symmetric. 



4 Surprises of the resolvent 

Consider a Fock vector $ e h = 7i <8> r 5 (£ 2 (-K)) belonging to the range of the 
resolvent 

$ = Rph <g> V>(«) = / dte"' rf C/' t /i(8)^(v) = {$ n (a;)}, 



/o 
$„(•) : iR n -> W, w = {wi,...,w n } 

with components (3.13): 

$„M = J™ dtf^{-tJi(G + t)-L*wJ*v(-T)dT}<i> ni t(w), 



i x 7 

for commuting operators L, W, and G as above: 

L= J dE x L(\), W = f dE x W(X), G= J dE x G(X), 

L(X) = 2tp(X)e~ l ' s>w (2 - iX)' 1 , W(X) = e lZ(X \ G(X) = -iu(X) + p(A) 2 /(2 - iX). 

Denote by 0„ ;t the Fourier transform of the function <f) ntt (uj) with respect to variables 
u = {ivi,...,u n } : 

4>n,t{r) = f[((W - I)I[ ,t)(T k )v(r k -t)+ v(r k - t) + LI m (r k )]h, 

where r = {ti, . . . , r n }. Let /C be a subset in {1, . . . , n} and let K a be its complement. 
Put 

p lc%) =U((W- I)v{r k -t) + L)l m (r k ) e B(h). 

then 

&m(t) = £(Vg(r) II «fo» " *)V (4-1) 



£ x m&K a 
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The functions PJct( T ) have discontinuity in hyperplanes where variables r k change 
the sign: 

lim P^(r) = I Ka (k)P^(r), 

>( n )rfV\ _ T /MpWu I I (\\T T\7?J +\ I T \ d("-1) 



T lim o ig'M = I^(k)P^>(r) + (^ - /)«(-«) + LjP^ >4 (r). (4.2) 

Therefore, from (4.2) follows: 

Ptt(r)\?=- + °o = -((W-IH-V+Ly^l^ljcik). (4.3) 

Let us find values of jumps of n ,i( r ) when r k changes the sign. Note that 
lim </y t (r) = v(-t)^ n -i it (ri, . . . ,T k -i,T k +i,. ■ ■ ,r n ). 

Taking into consideration equations (4.1) and (4.3) we obtain the amplitude and 
phase jumps for functions from the range of the resolvent of the unitary group Ut : 

lim 4> n ,t{ T ) = w lim 4>n,t{r) H-L0 n _i t (ri,...,T fc _i,r fe+ i,...,r n ). (4.4) 

Denote by T>w,l — D <S> T s {Wl(lR \ {0})) a vector subspace in h that satisfies the 
condition (4.4). By A(5±), A(<5±), and N we denote operators acting on symmetric 
Fock vectors as follows: 






K\{0})' 



x JJ rfr m ($ n , *„) w (ri, . . . , r fc _i, e, T fc , . . . , r n _i), 



A(«J±*) n (r) = lim £ #„+i(ti, . . . , r fc _i, e, r fc , . . . , r n ), N* n (u) = n* n (w). 

e ^ fc=i 

The boundary condition (4.4) in this notation looks as follows 

(iV + 1) _1 (7 <g> A(<J+) - W <g> A(<L))tf = (L (8) /) *. (4.5) 

Let us prove that the operator 

H = iG®I + I®V+iL*W® A(8-), E = F* [ dTa + (T)a(r) id T F^ r , (4.6) 

J]R\{0} 
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is symmetric in T>w,l- Let $, \P e 2}w,l an d let B be a Hermitian operator so that 
domB ®/D T^w,l- Integration by parts implies an identity, where the difference of 
substitutions is expressed through operators A(5±) : 



($,5<8>Etf)-(£<8>E$,tf) = if$, B ® (A(<L) - A(«J+))* J (4.7) 

Using boundary condition (4.4) for functions <3> n and 4f n we find the difference 
between the substitutions (4.6): 



ii $, B <g> (A(5+) - A((L))tf) = z($, (W*BW - B)A(<L)#) 

+i(L $, 5L *) + i(W A(<L) $, 5L *) + i(L $, 5W A(<L) #). (4.8) 

In the particular case B = I, the equation (4.8) becomes simpler: 

i($,(A(<J + )-A(<L))tf) = i($,L*L^)-(iLW®A((L)$,^) + (<l>,iLW®,4((L)^). 

Now the identity %G — iL*L = (iG)* and equation (4.7) prove that the operator H 
is symmetric on T> WL : 

($, K9) = ((/ ® E) $, *) + ($, {iG <g> / + iL*W <g> A(<L)}tt) 

-i($, J <g> (A(<J+) - A(<L))tt) 

= ((/ <g> E) $, #) + ($, (iG)* <g> /*) + (zL'W <g> A(<L)$, *)) = (H$, *). 

Let us find how the generator of the unitary group Ut acts on Fock vectors 
belonging to the range of the resolvent. Let f 6 h, $ = R^h <g> ip(v). From the 
definition of the generator we have: 

— 1 rf 1 r°° °° 1 r 

(*,H$) = lim T -($,^) = T / cftE-r / dr 



n=0 



x (&(r), ^e-( G ^-^-^^r^(-)^0„, +s ( r) ) 



w 



s=0 



(4.9) 



Note that the functions 4> n ,tij) depend on differences r k — t. Hence, 

d ~ n <9 

at ~ dr k 



d ~ n d ~ 

\t(r) = - £ ^-<j> n , t {r) = iF^E^r). (4.10) 
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On the other hand, from the definition of the operator A(5-) we have 

J^ T (A(<L)0 nii )(T) = nv(-t) <t> n -l,t{T). (4.11) 

Now from the (4.9) and form equations (4.10), (4.11) we obtain 

(*,H$) h =/ dt(rk,iGe-< a+ '*- iLW fo«- r )* r h) H + dt £ ^ A *„ <*t 
./o io ^ n! j(iR\{o}) 



x(Vn(r),e- (G+ ^^ L * H/ Jo^-)^( i G' + iL*Wv(-t) +2 V . 



() N v*(r; 



fc=l 



*, jiG + i^W" <g> 4(<L) + J <g> e}$) , 



that is the generator H of the group U t coincides with the generator (4.5). Thus we 
have proved the theorem. 

Theorem 4.1 The generator H = iG®I + I®E + iL*W ®A(5-) of one-parameter 
unitary group Ut is symmetric in T>w,l- 

We have proved that the operator H is symmetric without the assumption that 
operators L, G, and W commute. There is no conceptual difficulties to extend the 
construction of a symmetrical boundary value problem to a wider class of generators 

H = iGH-/<8E + ^L|WV®An(<5-) (4.12) 

with the boundary condition 

(N+ 1)~\I <g> A £ (6 + ) - J2 We, m ® A m (5-))^ = (L e <g> /) *, (4.13) 

m 

where W = {W £m } is (M x M) is a unitary matrices with coefficients from B(7i), 
{A e (g) : g e L 2 (M), 1 < t < M} are annihilation operators T 5 (L 2 (iR M )) which 
commute for different indices £, and G = — iHq + \Y^ U\L £ . 

5 The Markow evolution equation 

Let us consider how the Markow evolution equation can be derived from the 
boundary value problem for the Schrodinger equation 

ltf(t) = f-G + il ® E - ]T L *iWi, m ® A m (<L) W) 
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with the boundary condition (4.13). Let B be a Hermitian operator from B(7i) and 
let h, g G D. Consider an equation for mean values (g,P t (B)h)n = (U t g<S>^(0)\B<S> 
I\U t h <g> ^(0)) h . From (4.12) we have 

^-{g, P t (B)h) n = -((G + ]T L* t W e , m A m (S.))U t g ® *(0)|S ® 7|C/ t /i <g> tf (0)) 

-(£/, 5 <g> tt(0)|£ <g> /|(G + £ L|WV An(5_))C/ t fc ® *(0)) 

C/i (? <g) *(0)|S ® (A(<5+) - A(<L))|J7 t /i ® *(0))Y (4.14) 

Now the equation (4.8) reads as 

U, B ® ((A(S + ) - A(<L))tt) = £($, (^ m B% " B)A m (6-)v) 

+ Y,(Le®,BL^) + Y,((We, m A m (6_)$,BL e y) + (Lz®,BW e>m A m (6_)y) 
t e,m ^ y 

(4.15) 
Since A(<L)#(0) = 0, A(<L)tt(0) = 0, (4.14) and (4.15) imply: 



j t (g,P t (B)h) H 



4=0 



(g, C(B)h) = -(Gg, Bh) - (g, BGh) + ]T(L^, BL t 



Hence, we obtain an infinitesimal map £(•) of the Markov evolution equation in the 
standard Lindbladien form: 

j t P t {B) = C(P t (B)), C(B) = -G*B -BG+ l -Y. L \ BL t 
for G = -iH + ±J2L* e L e . 

6 Concluding remarks 

Among the most important problems which stay unsolved in this paper we can 
mention 

-a proof of srs-convergence U t J t * — ► u(0, t) for Hamiltonians with noncommut- 
ing coefficients H, R, and K; 

-an extension of a theory to the case E = / \uo\ 2 a + (uj)a{oj)duu or E = / (uu\ 2 + 
c 2 ) 1/2 a + (uj)a(uj)duj, uj e R; 

-a study of conditions necessary and sufficient for the symmetric boundary prob- 
lem described in §4 to be essentially self-adjoint; 
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-a generalization of quantum stochastic calculus to equations with nonadapted 
stochastic differentials extending the Ito multiplication table. 
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